Approach to equilibrium in adiabatically evolving potentials 
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For a potential function (in one-dimension) which evolves from a specified initial form, Vi(x), to 
on : a different Vf(x) asymptotically, we study the evolution, in an overdamped dynamics, of an initial 

probability density to its final equilibrium. There can be unexpected effects that can arise from the 
time dependence. We choose a time variation of the form V(x,t) — Vj{x) + (Vt — Vf)e~ xt . For a 
. Vf(x), which is double welled and a Vi(x) which is simple harmonic, we show that in particular, if the 

(£) • evolution is adiabatic, this results in a decrease in the Kramers' time characteristic of Vf(x). Thus 

the time dependence makes diffusion over a barrier more efficient. There can also be interesting 
resonance effects when Vi(x) and Vf(x) are two harmonic potentials displaced with respect to each 
other that arise from the coincidence of the intrinsic time-scale characterizing the potential variation 
and the Kramers' time. Both these features are illustrated through representative examples. 
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Ch 1 I. INTRODUCTION 

o 
o 



The pioneering work of Kramers [| on thermally activated barrier crossing has provided an understanding of the 
microscopic mechanism underlying the Arrhenius temperature dependence of crossing rates Q. Several variants of 
the basic problem have subsequently been studied in the literature. A fair amount of attention has recently been 
devoted to the study of more complex nonequilibrium systems. These include the cases of diffusion over a barrier in 
the presence of harmonic force HllHISHISHIIland diffusion over a fluctuating barrier [Til HI llH 111 ITl 111 Il7| . 
£S| The hallmark of the former situation is the phenomenon of stochastic resonance , where the signal-to-noise ratio of the 
04 ■ system response to an applied harmonic force displays a local maximum as a function of the diffusion constant or the 
temperature. In fluctuating barriers, the discovery |12| | that the mean first passage time has a minimum as a function 
of the correlation time characteri zing the fluctuation has prompted a wide variety of investigations. The problem 
of surmountin g p otential barriers |l7l ITsl Il9l | has gained importance in other fields of science such as evolutionary 
computations |2Ct l2l| and global optimization as well. 

In the present work we consider the situation of barrier crossing of a time-dependent potential which adiabatically 
evolves from Vi(x) at t = to the potential Vf(x) as t — ► oo. In such a situation there will be an eventual equilibrium 
• distribution given by 

8" P e9 ~exp-ME), (!) 



where e is the diffusion constant and the approach to this equilibrium will be governed by a characteristic time, 
which differs from the corresponding characteristic time for the stationary potential V(x) — V/(x). We note that the 
characteristic time is the same as Kramers' time where the potential Vf(x) is one with a barrier . 

Our main results are given in the next Section, where we derive the time-dependent probability distribution for a 
specific form of the time variation leading from the initial potential V* (x) to the final Vf (x) . Other forms of the time 
variation can be treated by a simple extension of the techniques outlined there. In Section III, two specific examples 
of evolving potentials are considered. When the time-scale of the perturbation matches the Kramers time for the 
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stationary potential Vf(x), there is a resonance which delays the onset of equilibrium. This case is treated within 
the time-dependent perturbation theoretic method of Dirac. The second case we study is one where Vi(x) has a 
single minimum while Vf(x) is bistable. By reducing to an effective two-state dynamics, we show that the Kramers 
time for the stationary potential Vf(x) is reduced. The paper concludes with a summary and discussion in Section 
IV. Our result also sheds some light on the global optimization scheme recently introduced by Hunjan, Sarkar and 
Ramaswamy (HSR) Ullp. 



II. THE TIME-DEPENDENT DISTRIBUTION 



For concreteness, we consider the time-dependent potential 



V(x,t) = V f (x) + [V i (x)-V f (x)]e 



(2) 



which evolves via homotopy from Vi(x) at t = to Vf(x) at t — > 00. The Fokker-Planck equation for the probability 
distribution P(x, t) is 



which, with the substitution 23, 24. 2 



dt dx\ dx } dx 2 ' 



(3) 



P(x, t) = 4>(x, t) exp - 



2e ' 



reduces to 



90 
~dt 



Ho^ + H^e 



-At, 



where (primes denoting differentiation with respect to x), 

2 2e f 2e 

and 

= Vi(x)-Vf(x). 



4e 



Hq satisfies the eigenvalue equation 



H ip n = -E n tp n (x) 



(4) 

(5) 

(6) 
(7) 

(8) 
(9) 



with nonpositive eigenvalues. By construction, the ground state has eigenvalue E = 0, with the eigenfunction 
ipo(x) = Aexp— Vf(x)/2e, A being a normalization constant. Denoting the space-independent part of Hi by Vo(t), 
the solution of Eq. (0 can be written as 



Application of the standard techniques of the Dirac time-dependent perturbation theory leads to 

Cn(t) = J2cm(t)(m\H 2 \n)e-( E ™- E ^ 



(10) 



(11) 
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where 

H 2 =e- xt (Hx-V ). (12) 

The perturbative expansion for the coefficients in Eq. (|10l) is in powers of H 2 

00 

c«W =^c„,(t) (13) 
j=o 

with 

c„o = (14) 

and for j > 1, 



(*) =^2c md ^ 1 (t){m\H 2 \n)e-^- E ^ t . (15) 



It can be seen immediately that the c„o's are constants determined by the state of the system at t — 0. Assuming 
that the system is in the equilibrium state of the potential Vi(x) at t = 0, namely 



P(x,0) = A cxp- — , (16) 
e 



this implies that 



0(o;,O)=P(a ; ,O)exp^=^oexp-^ (17) 
with normalization constant Aq. The constants c„o are now determined from the intial condition, as 

f°° Vi 
c n0 = A / dxifjn(x)exp~7r -■ (18) 

J-oc ^ e 



Substituting this in Eq. I|15|l . the complete solution to the problem can be obtained using Eqs. I|13fl and (|10fl . 

Yl 

2e ' 



Note that as t — > 00, (/)(x,t) in Eq. (fTT)f) tends to cooip (x) = cooexp — and therefore, 



P(x,t ^> 00) Aexp (19) 

the equilibrium distribution corresponding to Vf(x). 

Also note from Eqs. (|11I12I) that since H 2 has the time dependence exp— Xt, there will be a resonance when 
E m — E n = A, giving a secular growth of the first order term, c n i(t) cx t. This is analogus to case of the time- 
dependent perturbation theory |2q . |2?| in quantum mechanics. 

III. APPLICATIONS 
A. Case I 

Consider first a case where the initial and final potential have the same number of minima. Specifically, we take 
Vi = (x — a) 2 and Vf(x) = x 2 , namely harmonic potentials that are spatially displaced. This leads to 



AV = a 2 - 2ax, 



(20) 
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Hi 



which has the space-independent part 



giving 



A(a 2 — 2ax) 
2~e 



2ax 
e 



T . . Aa 2 are 



H 2 {x,t) = — (A + 2)e 



-A/, 



The leading term in the expansion, namely 



d 2 i x 2 

H = e— + 1 

ox e 



has the eigenvalue spectrum E n — 2n(n = 0, 1, 2, . . .) with eigenfunctions 



2"n! 



-1/2 



H n (x/y/e)exp- 



2e 



where TL n {y) are the usual Hermite polynomials. The time-dependent probability is therefore 



t) — 4>{x, t) exp(- 



x 2 (a 2 - 2aa?) g _ At . 



2e 



2c 



giving, for <p(x,t), the expansion (cf. Eq. Q1U[0 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



^(x,*) =^c„(t)exp(-2nt- ^(l-e- At ))^„(a 

n 

To first-order in the perturbation expansion, we find 



(27) 



c„(i = 0) 



2et 



e 2e 



[2 n r 



(28) 



At the lowest order of perturbation theory, only the c„o's, which are given by Eq. (|28|l matter. Straightforward 
algebra now shows that 



P(x,t) 



1 

7T6 



exp - 



x 2 + a 2 



2aa;e 



e e 

,-At\ „-2t 



ax(e At ) - e 
exp — > '- + 0(c n i), 



(29) 



Since all terms to first order have not been included in the perturbation, the equilibrium distribution is not properly 
normalized as t — > oo and has the extraneous factor exp(— a 2 /e). After computation of c„i from Eq. 1)15(1 we find, 
after taking the appropriate matrix elements and carrying out the integration over time, that 



5 



C n l (t) 



a I n + 1 



1-e 



-XI 



c„+i(0) 



A + 2 \ a I n 
A-2 jel 2 



1-e 



-(A-2)t 



C„-l(0) 



(30) 



The coefficients c n ±i(0) are known from Eq. I|28|l and after some amount of algebra we get, correct to first order in 
the perturbing 'Hamiltonian', 



P(x,t) 



2axe 



-XL 



Aax(e 



(A-2)e 



(31) 



The effect of time-dependence in the potential can be seen by contrasting the above result, Eq. (|31|l with the 
sudden limit, when the potential is instantaneously changed from Vi(x) to Vf(x). The initial probability distribution 



corresponding to Vi (x) is 



and this approaches the equilibrium distribution corresponding to Vf (x) as 



P(x,t) 



1 



2axe 



(32) 



The coefficient =22 of e~ 2t has the extra factor (1 + 2/ (A— 2)). Note that the time-dependent perturbation effectively 
keeps the system from attaining equilibrium by always managing to cause transition to neighbouring states. The 
approach to equilibrium depends on the value of the adiabaticity parameter, A, and there are three different regimes 
of interest. 

• When A > 2, the approach to equilibrium is governed by e~ 2t but the coefficient of this term is significantly 
increased. 

• If A < 2, the approach is controlled by e~ xt and in the long time limit 



P(x,t) 



(ire) 1 / 2 



1 + 



4ax/e \ xt 
2-A 



(33) 



• Finally, we have the extremely interesting situation of A ~ 2, in which case 



„, 1 x 2 / Aaxt \ . 

P( M ) = ^e--^l + — e-» + ...J (34) 

This is the resonance that we have discussed already, which shows up as the coefficient of the usual correction 
to P eq (x) diverging with time. 

This divergence of the coffecient of e~ 2t in Eq.(34) would eventually get transferred to the argument of the 
exponential function as is usual in such cases. This can be explicitly verified in this case, because an exact 
solution for harmonic potentials has been written down by Hanggi and Thomas [28j. The answer for P(x,t), 
adapting the work of Hanggi and Thomas to this situation is 



P(x,t) 



7re(l 



1/2 



exp 



(x - ae- 2t [l + 2(1 - e-( A -2)*)/(A - 2)])' 



') 



(35) 



If we expand the exponential in powers of 'a' all the three cases cited above are exactly reproduced. This shows 
that the method of quantum mechanical time dependent perturbation theory that we have adopted here in 
capable of yielding the correct results. 
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B. Case II 



We now turn to a situation where an initially single well structure, Vi = x 2 , crosses over to a double well structure, 

2 4 

Vf = — + as t — >oo. The approach to equilibrium in double well potential is governed by the Kramers time, the 
long time scale coming from the possibility of noise induced hopping . Following the procedure outlined in Eqs. (|4I5|) . 

we get 



rf 2 (x*-xf l + 3x 2 
H ° + ~Te — (36) 



Hi = |(l-x 2 )- 1^(3-^(1-^) 

_ ± a *(6-x 2 ) + l-x 2 (3-X , ) 3 {l-e-* t ). (37) 

The low-lying part of the eigenvalue spectrum of Hq is characterized by a set of close doublets with exponentially 
small separations, while the gap between two doublets is of 0(1). The ground state Eq = 0, while the first excited 
state is the ground state of the supersymmetric partner of 



Vf) 2 Yl 

4e 2 



(38) 



and is exponentially small pfll l30j| . The next excited state has eigenvalue approximately 2, and hence we can treat 
the dynamics of the low lying states as that of a two level system. Denoting the two states by </>o and <^i, with 
eigenvalues and 5, then 

4>{x, t) = co(t)to(z) + ci(t)e- a Vi(aO (39) 
The dynamics of cq and c\ is governed by 

c = (0o|#i|0o)e~ At coW + (MHi\<j>i)c 1 (t)e~ ix+ V t (40) 

C! = (0 1 |i/ 1 |0 o )e- (A - 5) *c o (t) + (0 1 |iy 1 [^ 1 )e- A *c 1 (t) (41) 

Since the perturbation Hi is even, ((f>i\H\\(f)o) — 0, decoupling cq and c\ Integrating Eq. (|41() and dropping terms 
like e~ 2Xt which are unimportant for t > A -1 , we find 

c 1 = Cl (0 1 |i? 1 |^ 1 )e- At (42) 

where 

% = |(1 - x 2 ) + lx 2 (3 - x 2 )(l - x 2 ) - A^( 6 _ ^) (43 ) 

The primary contribution to (</>i |0i) comes from the vicinity of x — 1 since </>i is an antisymmetric wave function 
strongly peaked near x ± 1. Since ~ — we have 

Cl (*) =J B^_ex P -^ (44) 
with B a constant of integration. From Eq. I)39|l we find, after a series of standard manipulations, that 
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cf>o(x) 
1 + f(x)e~ s ^ 



(45) 



In the above a is a measure of the strength of the ground state wave function at the origin, and where 



The inverse of 5 e ff gives the effective Kramers time for the system and is shorter than the scale for the time 
independent system. This speeding up is most effect in the adiabatic limit, namely for A <C 1. 



IV. DISCUSSION AND SUMMARY 



The above result is a simple analogue of the global optimisation principle on an evolving energy landscape. In 
this case one is interested in finding the minima of a multidimensional potential energy surface which constitutes the 
energy landscape in problems such as protein folding or finding the ground state configuration of atomic or molecular 
clusters. The observation of HSR that continuously and adiabatically varying potentials assist approach to 

the desired configuration at t — > oo by avoiding trapping in local minima. We have shown in a model system, a similar 
time dependence , the decrease of Kramers time makes escape from a trapping potential easier. 

The above demonstration of a reduced Kramers' time is for a one dimensional system. The extension to two 
dimensions is reasonably straightforward following the technique in ref. [30] when a well defined tunneling path exists 
between the two minima. Extensions to more general situations and to three dimensions is being investigated. We 
have also seen that this technique of dealing with time dependent perturbations can model the stochastic resonance 
as a kind of parametric resonance. With the emerging importance of stochastic resonance in biological systems fit! ] 
it is possible that yet another way of looking at stochastic resonance can yield new insights. 
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